Maximal entanglement from quantum random walks 
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The conditions under which entanglement becomes maximal are sought in the general one- 
dimensional quantum random walk with two walkers. Moreover, a one-dimensional shift operator 
for the two walkers is introduced and its performance in generating entanglement is analyzed as a 
function of several free parameters, some of them coming from the shift operator itself and some 
others from the coin operator. To simplify the investigation an averaged entanglement is defined. 
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I. INTRODUCTION 



Quantum random walks [1-4] (QRW) are a generaliza- 
tion of classical random walks. At each step of the latter 
the position of a classical particle (walker) is shifted ac- 
cording to the result of tossing a coin. Thus, the position 
of the walker is decided by following a probability dis- 
tribution over classically well-defined positions. In the 
QRW case the tossing of the coin is substituted by the 
action of a unitary operator (that shall be called coin op- 
erator) on a 2-state system (called coin state) and the 
motion of the walker by a unitary transformation of the 
position state by a shift operator. The state of the walker 
is in general a quantum superposition of several position 
eigenstates. 

Many aspects of such systems have been studied in the 
past [5-8] . In [9] the consequences of making one walker 
move in more than one dimension are analyzed. The in- 
terplay between decoherence and entanglement has been 
unveiled in [10] while in [11] two coins have been used. 
The so-called meeting problem has been studied in [12] 
for two walkers. We end this short review of recent results 
by recalling Ref. [13] where the evolution of the chirality 
of the coin has been investigated. 

QRW also attract interest in some specific research 
topics. For example, the efficiency of the energy transfer 
in photosynthesis can be raised at the level of about 99% 
if it is modelled by a QRW interacting with the thermal 
fluctuations of the environment [14]. Besides, quantum 
transport properties of electron systems and dielectric 
breakdown driven by strong electric fields can also be 
studied by utilizing QRW [15]. In ancilla-based quan- 
tum computation schemes [16] the ancilla plays the role 
of the coin and the measurement produces entanglement 
between the qubits. Moreover, experimental set-ups for 
QRW are nowadays available [17-19]. 

In the present paper we propose a quantum walk made 
of one spin i state (the coin) coupled to two particles (the 
walkers) whose positions can take any integer value on an 
infinite line. QRW with two walkers are a valid tool for 
obtaining highly entangled states [20,21]. Indeed, per- 
forming measurements on the 2-state coin system after n 



random walk steps may yield position states with strong 
entanglement among the two walkers. In particular, our 
QRW has been modelled in such a way that a further 
position measurement on one walker fixes the position of 
the other, thus rendering the correlation among the two 
walkers a maximum. 

The relevance of entanglement in QRW is crucial. 
While this property does not appear in classical walks, 
it becomes the characterizing ingredient of QRW. The 
interest of entanglement lies in the fact that it is of ut- 
most relevance for research in quantum computation pro- 
tocols and quantum information [22-26] (cryptography, 
communication, algorithms, etc.). We want to study the 
growth of entanglement among the two walkers and, in 
particular, to discover whether maximal entanglement 
may be achieved. We also pay attention to the quan- 
tum probability P to obtain such entangled states after 
the coin measurement. 

As a novel ingredient of the present paper, we will in- 
troduce a shift operator that to some extent generalizes 
the ones so far used in literature. Its form contains sev- 
eral free parameters which, added to the parameters of 
the coin operator, supply enough freedom to easily find 
conditions under which entanglement is maximally en- 
hanced. 

The plan of the paper is as follows. In section II the 
model for QRW is presented in detail. It includes a dis- 
cussion about the general type of shift operator. An an- 
alytic study of the QRW during the first steps in the 
walk is shown in section III. This approach, though ap- 
plicable to a very limited number of random walk steps, 
turns out to be extremely useful to interpret the numer- 
ical results derived in section IV. Indeed, in the latter 
section, thanks to a numerical treatment of the random 
walk evolution, we study the entanglement at an arbi- 
trarily large number of steps. In this section an averaged 
entanglement will be defined to simplify the search for 
the conditions of maximal entanglement. Some conclu- 
sive comments will be given in section V. 
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II. THE SPIN-PARTICLE QRW 

We begin by stating the precise definition of the model 
under study. The two walkers can roam on a line of 
discrete spatial positions. The accessible position eigen- 
states are for i e Z (i = is the origin) where the 
first item refers to the position of the first walker and the 
second one of the second walker. Thus the two walkers 
are supposed to stay together during the walk. The coin 
state is a spin i system. Its eigenstates are the spin com- 
ponent along the Z-axis. The complete quantum state 
is a sum of terms of the form \(f>) = \s) ® where the 
first factor refers to the spin and the second one indicates 
the position state. The <g> symbol emphasizes the tensor 
structure of the state space. The most general pure posi- 
tion state contemplated by our model is \ip) — J2 i Ci\i, i) 
with ^ \a\ 2 = 1. The eigenvectors of the Z-component 
of the spin are denoted by an arrow: | t) (\ I)) stands for 
the Z-component +| (—5). Its most general pure state 
is d t | t) + d±\ I) with \d^\ 2 + \d±\ 2 = 1. A similar model 
was introduced in [20]. 

The above model is hardly realizable in laboratory 
(for a more physical model with two walkers see for ex- 
ample [27] and references therein). Rather it must be 
viewed as the limiting case where entanglement is ex- 
pected to be maximal. Indeed, imagine a more general 
QRW model for which the position state of the two walk- 
ers after the spin measurement (performed after n QRW 
steps) is c ij\hj) € Z and £\ ^ \c i:j \ 2 = 1) where, 
as before, the first (second) item in the ket refers to the 
first (second) walker. This state displays large entangle- 
ment among the walkers if after a further measurement, 
this time on the position state of one of the two walkers, 
the position of the other walker gets fixed with maximum 
probability. Evidently, the limiting case of the above cir- 
cumstance occurs when j is a function of i, j = f(i) 
where / is a predetermined known function. As for / 
we have chosen the identity function. Moreover, since we 
have chosen such an idealized QRW model, the precise 
connection among the spin variable and the walkers will 
not be better specified. In the rest of the paper we will 
study the entropy that entangles the two particles in the 
state X)i c iK;*) after the coin measurement. 

Since we want to study the generation of entanglement, 
throughout the paper the factor corresponding to the po- 
sition of the walkers in the initial state |^)o of the QRW 
shall be free of entanglement, 



(rf t |t) + ^U))®|o,o), 



(1) 



where the above normalization condition regarding coef- 
ficients df,di is assumed and the subscript stands for 
initial state (0-th step of the QRW). Furthermore, since 
for any spin state |s) there exists some axis Z' such that 
I s ) = I t')i the state in (1) can be simplified to (dropping 
primes) 



>o = |t>®|0,0>, 



by adequately choosing the reference axes. Such a simpli- 
fication does not entail a loss of generality for the reason 
that the position state and the coin state in (1) are un- 
correlated. 

A single step in the random walk consists in the con- 
catenation of two operators, t/ s hift • U co in, iteratively ap- 
plied to \4>). Operator t/shift acts on position states and 
£4oin on spin states. 

Let us begin by introducing the shift operator [7 s hift • 
A quite general form is 

E4hift = a n I tXt I £ I* + «, * + «><*, *l 



£4? 1 



t)U \^2\i + n,i + n)(i,i\ 



l) (l \ ^2\i + n,i + n)(i, 



+ /Z a u 



(«) 



(3) 



with oS™} complex coefficients. They must satisfy 
several constraints derived from imposing unitarity, 
E^shift^shift = I: (I is the unit operator). After some 
algebra, the l.h.s. of this condition is 



£ [ I 4-)(t I \i + n,i + n)(i + m, i + m\ 



( 



(n) (m)* 



(n) (m)* 



) 



n, t + n\ 

(n)* (m) 



+ I t)U I \i + m,i + m)(i 

I (n)* (m) 

+ 1 1) (t I \i + n,i + n){i + m, i + m\ 

( (n) (m)* . (11) (m)*\ 

+ I L)(i I \i + n, i + n)(i + m, i + m\ 

I in) (m)* , in) (m}*\ 

x ( a H a U +«1tHt ) 

which must equal 

i=(|t)(tl + U)U l)£M( 



{1,1 



(4) 



(5) 



To fulfil the condition, the first two terms in (4) must 
be zero and note also that the second term is the h.c. 
of the first one. These considerations suggest to define 



the complex 2-vectors V{ 



(n) 



(a 



(«) An) 

TT ' Ti 



>) and V 2 (n) = 



(a\J,a\,') because in terms of them the above con- 



ditions become v/ n ^ • V^™ 1 ^* 



(n) _ y{m)* 



0, v{ n) ■ v; 



r(m)* 



for all n ^ m and V 1 



(n) _ y(n)* 



y(n) y(n)* = y(n) y(n)* = 1 for aU n Sin(;e in complex 

2-space there are only two independent orthonormal vec- 
tors, there is no much freedom to choose and . 
The most general solution reduces the sum over n and 
m to one single term, n = p, m = q for p, q fixed inte- 



(2) gcrs and = (a, (3), V 2 {q> = {-/3*, a*) (for complex 



(9) 
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numbers a, /3) and all the other y-vectors equal to zero. 
Moreover the two coefficients a, (3 satisfy |a| 2 + |/3| 2 = 1. 
Thus, the shift operator is 

C/shift - (o| t><t 1 + 01 t)U l) Y^\i + P ,i + p)(i,i\ 

i 

+ (cm>UI-/m>(tl)£> + <M + <z><M|, (6) 

i 

with p, q e Z fixed and a, (3 verifying the above normal- 
ization. 

For p = q no entanglement is generated by the QRW 
as every term in U s hm would move the couple of parti- 
cles to the same spatial position. Instead, entanglement 
appears as far as p ^ q, the specific values of p, q being 
immaterial since different pairs p, q are related by rescal- 
ings on the line of discrete positions where the walkers 
wander. During the present study we will stick to the 
values p — +1 and q = — 1 all the time. 

Without loss of generality, we can take a real in (6). 
Indeed, the spin eigenstates | t) > I 4-) can be redefined 
in such a way to absorb the phase of a. Calling £ 
such a phase, transformations | t) e^/ 2 | t) an d 



I) 



-iC/2 I 



4.) eliminate it from (6). There is no simi- 



lar procedure able to withdraw the phase of f3 too. Thus, 
the [/ s hift operator depends on only two variables: |a| 
(henceforth called just a > 0) and arg(/3). In literature 
the choice a = 1, /3 = is generally found. 

The coin operator [/ co ; n can be represented by a uni- 
tary 2x2 matrix. The most general such a matrix is an 
element of the U(2) group, 



-VI -pe 



V 1 - p e 

y/P e 



-2rq 



(7) 



where < 9, 77 < tt and < tp < 2ir are arbitrary phases 
and < p < 1. When tp = rj (7) becomes an element of 
the subgroup SU(2). Actually, the overall phase tp plays 
no role in the following discussion, so in practice we could 
limit ourselves to study SU(2) matrices to represent the 
most general coin operator. 

Particular instances of the unitary matrix [/coin are the 
Hadamard coin (p = 1/2, tp = and rj = 9 = n/2), 



TT 1 ( +l +1 



(8) 



the Kempe coin [2] (p = 1/2, ip = rj = and 9 = ir/2), 
1 f +1 +i 



U K = — 



V2 \ +i +1 

or the Z coin (p = 1, r\ = ir/2, tp = and 9 arbitrary), 



(9) 



U z = 



+1 
-1 



(10) 



At any moment the walk can be stopped and the value 
of the spin measured. If the measurement was performed 
after each step in the iteration then a classical random 
walk would result. Hence, our interest for QRW consists 
in relegating the measure after many steps because in this 
fashion clear entanglement is displayed. Let us consider 
a measurement performed only after n QRW steps giving 
spin s and a position pure state, result of the wavefunc- 
tion collapse, denoted by |i/>)* and equal to 



(11) 



If the spin measurement yields up (down) spin we will 
write \ip)n P (IV^n ™)- The evident entanglement created 
among the two walkers associated with \ip) s n shall be cal- 
culated by the von Neumann entropy 



c,| 2 log 2 I 1 2 , 



(12) 



where the subscript indicates that the measurement has 
been performed after the n-th step. Certainly E n > 
and maximal entanglement is attained when the values 
|cj| are all equal. So, if the sum in \^) s n = X^ c *lM) 
contains N terms, (£ n ) max = log 2 N. We shall usually 
calculate the normalized entanglement defined as 



E n 



[En) 



(13) 



This choice of units allows to easily understand when 
maximal entanglement has been attained. As the num- 
ber N of terms in \ip) s n in general depends on the step n, 
also the denominator in (13) depends on n. When N = 1 
both E n and £ n vanish, the former due to (12) and the 
latter by definition. 



III. ANALYTIC STUDY OF ENTANGLEMENT 
GENERATION 

The first few steps of a QRW can be followed analyti- 
cally. Starting from state (2), using the coin operator (7) 
(where the irrelevant global phase e lv has been skipped) 
and the shift operator (6) it is easy to see that entangle- 
ment can only arise after two steps. Indeed, the walkers 
plus coin state after the first step (indicated with a sub- 
script 1) is 



|$i - {ay/p- /VT^e"^))! t) ® |1,1} 
-{fi*y/p + aVT~pe-^)\ 4.) ® I - 1, -1), (14) 

and the wavefunctions after the measurement, \ip)^ p = 
|1, 1) and \ip)i° wn = I — 1,-1), whenever they exist, are 
clearly free of entanglement. Let us assume that the spin 
measurement after the second step yielded +5. The re- 
sulting entangled position state is 



3 



|V>£ P cx (a^p - p^T~p e-^) 2 |2, 2) 



ay/l- p + P^fp e 



|0,0). (15) 



The presence of the proportionality symbol reminds the 
lack of a normalization constant. We will often write the 
wavefunction resulting from the spin measurement in this 
way purposely because it renders more transparent the 
evaluation of the probability P with which such a state 
comes out (P is just the square of the missing normal- 
ization). It is easy to verify that maximal entanglement 
is achieved when the condition 



P e 



ay/1- p + /Vpe 1 



(0+v) 



(16) 



holds. If the Hadamard coin were used, the condition 
would became \a + p\ — \a — P\ which is fulfilled when- 
ever a is real and P pure imaginary. For the Kcmpe coin, 
condition (16) reduces to \a + i/3| = \a — i/3| which is sat- 
isfied for any a and P as far as both are real. For the Z 
coin, condition (16) requires |a| = \P\. 

If the result of the measurement after two steps yields 
spin — i, then the position state is 

|V>>2° wn oc (ap{l - p)e- 2i ( e+ ") + \p\ 2 ^ p{l - p)e-' l{B+ ^ 
-a 2 x /p(l-p)e- i ^ - ap*p)\0,0) 
+ (ap*(l - p)e- 2 '" + (^)V/)(l-p) e i(M 



Vp(l-p)e- i(8+ "V 2i " 



— ap pe 



2i ")|-2,-2). 



(17) 



To find the circumstances under which entanglement be- 
comes a maximum we discuss different values of p sep- 
arately. When either p = or 1 then only one term 
survives in each coefficient of (17) and maximal entan- 
glement occurs for all 9, i], a and arg(/3). For p = \ 
entanglement is a maximum again for all possible values 
that the parameters can take. To see this we substitute 
p = \ in (17) and obtain after some algebra 

wn oc e-W+rt (\P\ 2 - a 2 + 2\a\p\ sin a) |0, 0) 

+e -2ir )e -iarg(« _ 1)cogA 

-i sin A) I -2,-2), (18) 

where A = arg(/?) — 6 — rj. It is easy to verify that both 
coefficients have the same modulus. 

An instance of maximum entanglement can also be de- 
duced for the general case < p < 1. First note that 
every term in the two coefficients in (17) are equal, apart 
from phases. Imposing that the phases are also equal 



tcrmwise (except for a possible global phase) , we discover 
that maximal entanglement occurs when arg(/3) = 8 + r/. 

A generalization of these analytic results to n iterations 
of the QRW is too complicated as the length of the coef- 
ficients in |i/))" p ' down grows exponentially with n. Then, 
in the following section we study the above problem nu- 
merically with a computer program. The program yields 
exact numerical results but without the corresponding 
explicit analytical expression. 



IV. ENTANGLEMENT AFTER AN ARBITRARY 
NUMBER OF STEPS 

In this section the entanglement created by measuring 
the spin of the coin after an arbitrary number of QRW 
steps having started with the initial state (2) shall be 
studied. To this end we introduce an averaged entangle- 
ment E n defined over the first n steps as follows. Consider 
n replicas of the QRW and evolve each of them indepen- 
dently of the others. In the first replica we measure the 
spin and then the entanglement after the first step ob- 
taining £i , in the second one after the second step obtain- 
ing £2, etc., thus gathering the collection £1, £2, £3, ■ ■ ■, 
£ n . On the understanding that each measurement yields 
the same result (i.e., cither all measurement outputs give 
up spin or all give down spin), we define 



1 " 
■ - V5 



(19) 



The first step (a = 1) is excluded from the sum in the av- 
erage because it generates zero entanglement in any case. 
Note that in (19) the subscript n does not mean the step 
after which the measurement is done, as in (12), but the 
number of steps over which the average is performed. 
The convenience of using such an averaged entanglement 
becomes obvious when it is 1 because, being the mean of 
many positive numbers not larger than 1, such a result for 
£ n necessarily implies that all of them are precisely 1, i.e. 
it identifies situations where entanglement is constantly 
maximal (therefore the precise value of n used to evalu- 
ate £ n is largely immaterial, as long as it is sufficiently 
large). This assertion applies equally well when £ n = 0, 
but clearly this case is physically less interesting. 

A computer code was prepared to implement the 
above-described procedure for the QRW during a num- 
ber n of steps. It calculates step by step the entangle- 
ment, either averaged or not, as well as the probability of 
having up or down spin after the measurement of the coin 
state and also the number of terms N in the collapsed 
position wavefunction (the latter enables us to compute 
the exact value of (£ B ) mas ). The numerical procedure 
used to extract N was demanding that \ci\ in (11) be 
larger than a pre-fixed threshold, in our case 10~ 10 . All 
the results presented in this section have been obtained 
by using this code. 
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The next three subsections contain the study for the 
three coin operators of section II, namely the Hadamard 
Uh, the Kempe Uk and Uz, respectively. The last sub- 
section is dedicated to the general coin case, (expres- 
sion (7) with tp = 0). 
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FIG. 1. Averaged entanglement obtained during the first 
200 steps of a QRW with the Hadamard coin (8) and shift op- 
erator (6) with a, P real. The plots for down (up) spin output 
of the measurements are shown with a continuous (dashed) 
line. Symbols | 1 /j) u P> down here and in the next figures refer to 
the output of the spin measurement. 



A. Hadamard coin 

This is the most frequently used coin operator for the 
kind of problems at hand. Its expression is given in (8). 
To study the entanglement generated by a QRW evolved 
with this coin, we first calculate £ n for n = 20 steps 
for real and positive values of a and f3 = +Vl — a 2 
in (6). The result is displayed in Fig. 1 as a function 
of a (0 < a < 1). 

The plots in Fig. 1 show that as a and (3 tend to be alike 
and real (this means both approach the value -j=) entan- 
glement tends to be maximal (minimal) if the measure- 
ment yields down (up) spin. The number of entangled 
terms N in the collapsed wavefunction after the measure- 
ment is equal to the step number n. This must be inter- 
preted as that entanglement gets richer and richer as the 
QRW goes on. However at the exact values a = (3 = 
entanglement suddenly drops to zero for both kinds of 
measurement outputs. This is how the vertical line in 
Fig. 1 (which applies to the down spin measurement) 
must be deciphered. Moreover, the probability of having 
such maximally entangled states (for a w f3 and both 
real) after the measurement is extremely small, P <C 1. 

The origin of all the above facts can be clarified by 
using the analytical expressions obtained in section III. 
First of all we recall that for any values of a and (3 the 
first step always gives no entanglement at all. For in- 
stance, taking a = f3 = -4= the state after one step of 



QRW is \c/))i = | t) <8> 1 1 , 1). Proceeding with the choice 
a = j3 = -j^, a calculation similar to that exhibited in 

section III allows to show that \<f>}2 = | t) ® 1 2, 2). By in- 
duction this can be immediately generalized to any step, 
getting \4>) n = | f)<g)|n, n). Hence, entanglement vanishes 
when a = /3 = 

Now we show why for step n = 2 the entanglement 
of the state resulting from a down spin measurement 
approaches its maximal value but suddenly drops to 
zero at a = f3 = Taking a, (3 in (17) real and 

close to each other leads to the expression |V , )2° wn 
(/3 2 - a 2 )(\ - 2, -2) — |0, 0)) which indeed displays maxi- 
mal entanglement although it suddenly becomes zero as 
soon as a exactly equals /3. This expression also indi- 
cates that despite the entanglement of the state resulting 
from a down spin measurement converges to the maxi- 
mal value, as a tends to {3 the probability of obtaining 
such a measurement output becomes very small, of or- 
der 0(\(3 2 — a 2 \ 2 ). It seems reasonable to expect that an 
analogous mechanism explains the low probability also 
for n > 2. 

In Fig. 2 the ratio £ n (this time non-averaged) for three 
different values of a (taking both a and {3 real and posi- 
tive) are shown as a function of n. It is clear that, as a 
approaches /3, both tending to the entanglement gets 
stuck to the maximum for more and more steps. There- 
fore we infer from Fig. 2 that when a = 0.71 the averaged 
entanglement after 200 steps is £200 ~ 1 while after more 
steps, say 800, it is patently less than 1, £ 8 oo < 1- In- 
stead for a = 0.7071, a figure closer to both £ 2 oo 

and £ 80 o are firmly anchored to 1. Recall that, although 
not shown in this graph, we have proved before that at 
exactly a — (3 = the entanglement drops to zero. 
The plot for other values of a are qualitatively analogous 
to the one shown in the jagged line of Fig. 2 (here and 
throughout the rest of the paper, the number 0.37 will 
just represent an arbitrary value of a). Thus, a method 
to obtain highly entangled states would be using real val- 
ues of a and /3, very close to each other, and look for a 
down spin measurement output in the first QRW steps. 
The difficulty is that, as explained above, such down spin 
outputs are very unlikely. 

Next we add a phase to f3 while keeping a real. In 
this way we will have covered all relevant (real or com- 
plex) numerical values for a, (3. No dramatic changes in 
the entanglement generation are found. In particular, no 
cases with arg(/3) 7^ exist with maximal entanglement, 
see Fig. 3 and Fig. 4. When the dependence on 

arg(/3) is qualitatively as in Fig. 3. For a — — (3 = 
the averaged entanglement goes to zero, as Fig. 4 shows 
manifestly, meaning that the true, non-averaged entan- 
glement vanishes at all the QRW steps. 

The use of the averaged entanglement £ n after n QRW 
steps must not divert our attention away from the fact 
that the real interesting quantity is the normalized en- 
tanglement £ n achieved by the quantum state after each 
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single measurement. Our computer code allows to dis- 
cover isolated cases of maximal entanglement (being iso- 
lated, they get lost in the average of £ n ). They occur 
mainly in the first few steps of the QRW and are quite 
likely since the probability of obtaining the related spin 
measurement is not very low. 

Let us enumerate all cases of this kind found when us- 
ing the Hadamard coin. We start by the position states 
originated by measurements that led to down spin: the 
entanglement is a maximum at the second step for all a 
or (3 (real or complex) (except for a — (3 — as we 

already know) and at steps n = 3 and 4 when the argu- 
ment of f3 is arg(/3) = ir/2, 3tt/2. If instead the output of 
the measurement at n — 2 is up spin, then the entangle- 
ment is a maximum for all a real and arg(/3) = tt/2, 3tt/2. 
When n = 2 these features can be proved by resorting 
to expressions (15) and (17) (but the formalism of sec- 
tion III does not suffice to investigate the presence or ab- 
sence of maximal entanglement beyond the second QRW 
step). Moreover, in all cases the probability for obtaining 
the indicated result from the spin measurement is around 
P <~ 0.2 — 0.5 and the number N of entangled terms in 
the wavcfunction collapsed after the measurement coin- 
cides with the QRW step number, N = n except when 
n = 3 for which N is 2. There seems to be no other steps 
n with £ n equal to 1. We have checked this assertion up 
to n as large as 1000. 



B. Kempe coin 



When using the Kempe coin (9), the average entan- 
glement £ n turns out to be independent of a, [3 as far as 
both coefficients are real and positive (but depends on n 
and on the result, up or down spin, of the measurement). 
However, the dependence on the argument of a complex 
j3 displays the rich structure shown in Fig. 5. The first 
aspect of this graph to highlight is that, once fixed the 
result of the spin measurement, the averaged entangle- 
ments at real and positive a, (3 (the rightmost or leftmost 
edges in the figure) arc indeed the same for the two values 
of a that have been utilized. The second aspect is that 
for down spin measurements and a = \(3\ = there is 

maximal entanglement as arg(/3) tends to become 37r/2 
but at this precise value, it suddenly drops to zero (indi- 
cated by the vertical line at this value of arg(/3)), a phe- 
nomenon similar to the one described for the Hadamard 
coin in Fig. 1. In fact, at a = ^ and /3 = — -±= the 

wavefunction after n steps is \cj)) n = | t) ® \n,n). As 
also happened for real values a ~ j3 in the Hadamard 
coin, the probability of having this output after the spin 
measurement is negligibly small for arg(/?) s=y 37r/2. 
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FIG. 2. Entanglement obtained in the first 800 QRW steps 
after a down spin measurement with the Hadamard opera- 
tor and a, P real. The continuous line shows the result for a 
value of a very close to = 0.7071067812 • • •. The dashed 
line corresponds to a value of a not so close to and the 
lower continuous (jagged) line to a value decidedly different 
from ^. 
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FIG. 3. Averaged entanglement £200 as a function of the 
phase of parameter (3 for a value of a far from with the 
Hadamard coin. 
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FIG. 4. Averaged entanglement £200 as a function of the 
phase of parameter /3 for a — with the Hadamard coin. 

The behavior near arg(/?) ~ n/2 can be described as 
follows: the averaged entanglement becomes £200 = 0.5 



G 



for a = \/3\ = -j= and arg(/3) close to 7r/2 while it drops 
to zero when the phase attains this precise value. This 
pattern can be spelled out by studying the non-averaged 
entanglement: for a s=s s=s and arg(/3) = 7r/2 it 
is maximal for even QRW steps and a down spin mea- 
surement while is zero for odd steps. Instead, for up 
spin output, the order is reverted: for even steps is 
zero and for odd steps is maximal. This explains the 
value 0.5 that £200 attains for both spin measurements 
near arg(/3) = 7r/2 in Fig. 5 (in other places where £200 
equals 0.5 the above-described altcrnance between even 
and odd steps is not seen). Unfortunately the probabil- 
ity of having such measurement outputs (down spin for 
even steps or up spin for odd steps) is negligibly small, 
P « 1. Furthermore, in all these cases the number of 
terms in the collapsed wavefunction is just N = 2 for 
any step n denoting a rather poor entanglement. When 
exactly a = \/3\ = and arg(/3) = 7r/2, the entan- 
glement vanishes at all steps for any measurement out- 
put (actually, the wavefunction after n QRW steps is 
\4>)n — I 4-) ® I — n i ~ n ))- This accounts for the zero of 
£200 in Fig. 5 at strictly arg(/?) = 7r/2 and a = |/ 



1 
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Again we study the non-averaged entanglement during 
the initial steps fo the QRW. If the measurement output 
is down spin then entanglement is maximal at the steps 
n = 2,3,4 of the QRW for all a, (3 real and at the step 
n = 2 for a real and any complex ft (the cases for n = 2 
follow also from (17)). When instead the measurement 
gives up spin, the QRW produces a state with maximal 
entanglement only at the second step for all a and /3 real 
(it matches with the discussion after expression (16)). 
All the above events with maximal entanglement have 
a significant probability to occur, P <~ 0.2 — 0.5 and 
the number of terms N in the collapsed position wave- 
function coincides with the QRW step n except for n = 3, 
down spin and a, (3 real, for which N is 2. Up to n = 1000 
we have found no other maximum non-averaged entan- 
glement results for any coin measurement output. 



C. Z coin 



In Fig. 6 we show the averaged entanglement obtained 
with the use of the Z coin (10) as a function of a for real 
j3. From this plot it is apparent that there are no val- 
ues of these parameters that allow to have maximal or 
almost maximal entanglement during the entire QRW. 
Moreover, having fixed the (real) values of a and |/3|, the 
plot of £200 as a function of arg(/3) (not shown) comes 
out flat, indicating no dependence on this phase. There- 
fore we conclude that the averaged entanglement never 
attains the maximum when the Z coin operator is used, 
implying that there are no values of a,/3 that allow max- 
imal entanglement for all (or almost all) n > 2. 
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FIG. 5. Averaged entanglement £200 as a function of the 
phase of parameter ft for various values of a with the Kempe 
coin. The thick continuous line refers to down spin measure- 
ment output with a = \f}\ = the thick dashed line to up 
spin and a = |/3| = the thin continuous line to down spin 

with a = 0.37 and |/3| = — a 2 ; and the thin dashed line 
to up spin with a = 0.37 and \/3\ = — a 2 . 
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FIG. 6. Averaged entanglement £200 as a function of a for 
real a.B with the Z coin. 



Maximal entanglement can be achieved only in spe- 
cific cases. A detailed analysis performed with the non- 
averaged entanglement enables us to discover that in- 
deed this event happens for all real a and all complex 
j3 when the measurement output yields down spin at 
the step n = 2 of the QRW. The number of entan- 
gled terms in the collapsed wavefunction |V')2 0wn after 
the spin measurement is N = 2 and the probability is 
P ~ 0.2 — 0.5, again definitely far from zero. Further- 
more, for a = = 4= and arbitrary arg(/3) maximal 
entanglement appears for n = 2,3,4 steps and with a 
similarly high probability. In this case, though, the num- 
ber N of entangled terms in the resulting wavefunction 
after the spin measurement is equal to the step n only 
for n = 2, 4 while for n = 3, N is 2. On the other hand, 
if the measurement output is up spin, then the entan- 
glement will become maximal only when a = \/3\ = 

and arbitrary arg(/3) at the second step with P = 0.5 
and N — 2 entangled terms. As a and \(3\ tend to 
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be equal, the entanglement grows and the probability 
approaches 0.5. All these results for n — 2 can be 
easily recovered by using the formalism of section III. 
In particular, we have |V>)2 P oc a 2 |2,2) + |/3| 2 |0,0) and 
|V>)2° wn oc a*(3*\ - 2,-2) - a/3* 1 0,0). 



D. The general coin 

We have also run a battery of QRW with the general 
coin operator (expression (7) with ip = 0) in order to 
look for precise combinations of the parameters in [/coin 
allowing maximal entanglement. The initial state was 
still given by (2). 

Two kinds of runs were done. In the first one the aver- 
aged entanglement was calculated after 200 QRW steps 
for various values of p, 9, rj (in [/coin) and a, arg(/3) (in 
£4hift)- The ranges being < p, a < 1, < 9, 77 < 7r and 
< arg(/3) < 27r, data were taken by varying the values 
of each parameter by jumps of 0.05. The goal was to 
find QRW presenting maximal entanglement during all 
the steps. 

The second kind of run was designed to find isolated 
cases of maximal entanglement during the first few steps 
of the QRW. For this reason the QRW were evolved only 
for 10 steps and the non-averaged normalized entangle- 
ment was extracted. The same five parameters as before 
were swept within their ranges. 

No highly probable parameter regions producing maxi- 
mal or near maximal averaged entanglement (£200 > 0.99 
and P > 0.15) were found. Instead our computer runs 
were plenty of isolated cases with maximal entanglement 
at steps n = 2, 3, 4 of the QRW, (S n = 1 and P > 0.15). 
There were so many of them that giving a complete 
record lies beyond the limits of this short paper. Most 
of those at step n = 2 are summarized in the general 
discussion for arbitrary p after (17). 



V. CONCLUSIONS 

A quantum random walk (QRW) model with two walk- 
ers has been devised to study the generation of maximum 
entanglement among the two walkers by the process of 
coin measurement. The model is rather unphysical but it 
has been formulated in the way it is in order to maximize 
the entanglement among walkers. 

The allowed positions of the walkers are the set of (pos- 
itive, negative or null) integers along an infinite line. The 
2-state coin is represented by a spin \ system. Every 
QRW was started with the state | f) ® |0, 0), that is, the 
spin in the (Z-component) up eigenstate and the walk- 
ers at the origin of the line. The only consideration used 
to select this choice of the initial state was starting the 
walk with a position state free of entanglement. With 
this premise, and by a redefinition of the origin, any ini- 
tial position eigenstate is as valid as any other. The spin 



state is also completely general because by an adequate 
definition of the axes orientation any spin state can be 
viewed as the up spin eigenstate along the Z-axis. 

We have introduced a shift operator in the QRW with 
two walkers. It includes two free parameters that, to- 
gether with those present in the coin operator, provides 
enough freedom with which a measurement performed 
on the coin state may yield maximal entanglement on 
the resulting position quantum state. We have studied 
the problem both analytically along the first steps of the 
QRW and by a numerical computer code that allows to 
probe an arbitrarily large number of QRW steps. Besides 
the entanglement, we have measured the probability of 
having such highly entangled states and the quality of 
their entanglement, given by the number of terms that 
the collapsed position wavefunction contains. We have 
also devised an averaged measure of the entanglement 
that simplifies the search for a QRW with all or almost 
all steps presenting high entanglement. 

We have used the Hadamard, Kcmpe and Z coin oper- 
ators as they are the most generally seen throughout the 
literature. Moreover, in subsection IV D we also tackled 
the general coin operator. 

In general maximal entanglement is more frequently 
generated when the spin measurement output is oppo- 
site to the one used as initial state, in our case down 
spin. 

We have found two kinds of situations with maximal 
entanglement: (i) near maximal entanglement through- 
out almost the entire QRW, regardless of the number of 
steps, but with negligible probability to occur P <C 1 
(the larger is the number of steps, the lower is the value 
of P) and (ii) exactly maximal entanglement during a 
few steps at the beginning of the QRW with moderately 
high probabilities, P ~ 0.2 — 0.5. 

There are no highly probable entangled states contain- 
ing a large number of entangled terms. Only during the 
first steps of the QRW one can achieve maximal entangle- 
ment while having a reasonably high probability. How- 
ever, such states contain a limited number of entangled 
terms precisely because they occur at the first steps of 
the QRW. This happens for instance for down spin mea- 
surement at the second step for arbitrary a and arg(/3) 
with any of the three coins mentioned in section II. The 
largest absolute entanglement with an acceptably high 
probability and with the above three coins is achieved 
after a down spin measurement, occurring with probabil- 
ity P = 0.25, at the n = 4 step and managing to entangle 
N = 4 terms in the position wavefunction (thus giving 
E4 = log 2 4 = 2). They arc: (i) Hadamard coin for all 
a and arg(/3) = 7r/2, 37r/2, (ii) Kcmpe coin with any a, ft 
real and (Hi) Z coin with a = |/3| = -4= and all arg(/3). 
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